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ABSTRACT 

The  thermal  theory  of  laminar  flame  propagation  for  hyd3:%>gen- 
bromiine  mixtures  is  described.     The  method  of  analysis  follows  the 
earlier  work  of  von  Karman  and  Mill^n  and  of  von  Karman  and  Pen- 
ner.     The  problem  is  materially  simplified  by  introducing  the  steady- 
state  approidjnation  of  classical  chemical  kinetics  for  tlie  concentra- 
tions of  H-  and  Br -atoms. 

The  general  formulation  of  the  problexn  is  presented  in  Sec- 
tion II. 

Approximate  solutions  of  the  relevant  iXiathematical  problems, 
utilizing  various  procedures  developed  by  von  Karman  and  his  collabo- 
rators, are  given  in  Sections  III  to  V  for  bromine-rich,  stoichiometric, 
and  hydrogen- rich  znixtures,  respectively. 

Numerical  values  for  the  physico-chemical  parajrneters,  which 
are  required  for  the  calculation  of  absolute  values  for  the  burning 
velocities,  are  considered  in  Section  VI. 

Thj|  calculated  burning  velocities  are  compared  in  Section  VII 
with  experimental  results  obtained  by  R.  C.  Anderson  and  his  col- 
laborators.    Reference  to  Section  VII  shows  that  the  variation  of  the 
observed  burning  velocities  with  mixture  ratio  is  predicted  correctly 
by  the  theoretical  calculations.    Absolute  values  for  the  laminar 
burning  velocity  cannot  be  estimated  with  certainty  because  of  the 
possible  existence  of  large  errors  resulting  from  extrapolation  of 
low -temperature  kinetics  and  heat  conductivity  data. 
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I.    INTRODUCTION 

In  rect^nt  years  a  large  number  of  publications  have  appeared 
on  experimental  measiirernents  and  theoretical  studies  of  laminar  flame 

propagation  in  prennixed  gases.    For  a  complete  bibliography  reference 

(J)  (?) 

should  be  oiade  to  the  books  of  Lewis  and  von  Slbe^  '  and  of  Jost^"'  or 

(3) 
to  a  survey  paper  by  Evans^  '.    ^Experimental  nneasuremcnts  on  hydrogen- 

(4   5) 
bromine  fla ones  have  been  described  by  Anderson  and  his  collaborators^    *    , 

As  yet  incomplete  theoretical  studies  on  this  flame  have  been  published  by 
J.CoHirschfelder  and  his  collaborators^    '    ',  who  are  attempting  to  carry 
out  a  rigorous  numerical  solution  of  the  complete  flame  equations  for  the 
stoichiometric  naixture  ratio. 

The  method  of  approach  which  is  followed  in  the  present  investiga- 
tion differs  in  scope  and  purpose  from  the  type  of  calculations  carried  out 
by  KirschfeH^r .    We  follow  the  work  of  von  Karman  and  Xlillan^  '  and  of 
von  Karman  and  Penner^    *        '  in  that  we  shall  content  ourselves  with  an 
approximate .:5olution  of  the  flanne  equations  for  a  thermal  theory,  on  the 
grounds  that  ^e  required  physico-chemical  parameters  are  not  known 
with  sufficient  accuracy  to  justify  numerical  solutions.     On  the  other 
hand,  in  order  to  obtain  a  reasonably  severe  empirical  check  of  the 
physical  model  and  of  the  mathematical  approximations,  complete  calcula- 
tions as  a  function  of  mixture  ratio  are  carried  out.    In  order  to  obtain  the 
desired  results,  we  shall  utilize  the  steady-state  approximation  of  classical 

reaction  kinetics  for  the  chain  carriers,  as  suggested  for  flames  by 

z'      /  (9) 

von  Karman  and  Penner     ';  although  the  steady-state  approximation  cannot 

be  justified  for  Br  atoms  at  elevated  temperatures,  it  appears  that  the 
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calculated  burning  velocities  are  nevertheless  approximately  correct*  . 

There  is  a  well-known  cold  boundary- value  difficulty  whicl 
enters  into  the  solution  of  the  one-dimensional  laminar  flame  equations. 
"We  shall  not  discuss  this  problem  except  to  note  that  a  convenient  rrsethod 
for  avoiding  complications  is  to  utilize  the  concept  of  an  ignition  temper- 
ature^    '    . 


*  This  result  can  be  established  by  an  extended  analysis  camparing 
reaction  rates  with  and  without  the  steady-state  approximation.     The 
relevant  analysis  has  been  carried  out  by  von  K^rm^Tn  and  Millin, 
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II.    GENERAL  EOUATICNS  FOR  THE  THERMAL  THEORY  OF 

LAMINAR  FLAME  PROPAGATION  IN  HYDROGEN- BROMINE 
FLAMES  WITH  THE  STEADY-STATE  APPROXIMATION  FOR 
H    AND  Br 

For  the  theoretical  study  of  laminar  flame  propagation  in  hydrogen- 
bronnine  flatxies ,  we  shall  utilize  the  accepted  reaction  scheiriC  for  isother- 
mal production  of  HBr  in  stationary  systems  at  temperatures  of  the  order 
of  500   K.    The  required  equations  ar^r     ' 

X'  +  Br2   ilU      2Br    +    X'     ,  (1) 

Br  +  H^      ^2^      HBr    +    H      ,  (2) 

H    +    Br^  JX     ^^^    +  ^^     »  (3) 

H     +  HBr^      H^    +  Br       ,  (4) 


2Br  +  X'    ^      Br-  +  X'       .  (5) 


Here  X*  denotftg  the  "third  body",  i.e.    any  of  the  chemical  species 
H,  Br,  *ij,  Br,,  or  HBr.     The  specific  reaction  rate  constants  k.  to 
kc  for  the  five  reaction  steps  are  indicated  above  each  of  the  arrows. 

It  is  well-known  that  the  steady-state  approximation  for  H-  and 

(12) 
Br-atoms  leads  to  the  relation^     ' 

k     1/2  1/2 

2k2(r-i)        (H,)  (Br   ) 

^mH^.  ^     k         - (6) 

^dt  ^4    (HBr)  ' 
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where  the  chcfTjical  symbols  in  parentheses  denote  concentrations  for 
the  chemical  species  and  t    represents  the  tii-ne.    For  thermal  flame 
theories  it  is  convenient  to  introduce  mole  fractions,  X;  we  identify 
HBr ,  Br-,  and  Hj  by  the  subscripts  1,2,  and  3,  respectively.     Then 

X, 


(H2)=j5^X3=|j(X3_^-^)  (7) 


and 


where  p  denotes  pressure,  R  is  the  molar  gas  constant,  T  is  the  temp- 
erature, and  the  subscript  o  identifies  initial  conditions. 

For  flames  without  inert  diluents  it  is  convenient  to  identify  the 

initial  composition  by  the  parameter  a.    let  X-       =  a,  X_        =  1-a,     Thei> 

"^  J  » o  ^  ,  o 

Eq.  (6)  becomes 

,      1/2         3/2  „  Y    3/2 

^         d(XJ       2k2t)       (^)       (a  . -j-Hl-a--^) 

RT    -at^ ;    IT  •  ^ ' 

(l-a)-(i-^)X 
If  ¥/.  ,  ^2'  *^*^  ^ X  denote  the  molecular  weights  of  HBr,  Br^* 


and  K,!  respectively,  then 

d(X,) 


W, 


1'    p 


1""3F"  RT 

represents  the  rate  of  change  of  the  mass  of  bromine  per  unit  volume 

with  time.    But  this  quantity  is  also  equal  to 

WjdXj 

if  m   is  the  nnass  flow  rate  per  unit  area  of  the  tr.ixture  (i,e,  the  eigen- 
value for  the  problem  of  laminar  flame  propagation) ,  x  equals  distance , 
and 
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-^  =Z    WjXj  (10) 

is  the  average  aiolecular  weight  of  the  gas  mixture.    Hence  Eq«  (9) 
becomes  k     ^  ^^  ^^^         X  X    ^^^ 

J2._JL  = ! , — .  (11) 

In  the  absence  of  diffusion  the  general  energy  equation  for  one- 

(9) 
dimensional  flow  can  be  written  in  the  form     ' 

4_    f  =  (9  -  1)  +     ^'yj-^'"!  {12) 

where    A    equals  the  thermal  conductivity j    the  reduced  tetrperature    is 
0  =  T/T,,  with  the  subscript  f  identifying  equilibrium  conditions  at  the 

adiabatic  flame  temperature;  Y ,  is  the  weight  fraction  of  the  j'th  chemical 

o 
species;  h  .  equals  the  standard  specific  enthalpy  of  formation  at  ternper- 

ature  T     of  the  j'th  species;  the  average  specific  heat  c     is  giveu  oy  the 

o  p 

•      C9)> 
expression  ^  '    ?« 


=  -_-J-_    ^   (Y.   .  -  Y.      )h°   .  (13) 

P  (I-Oq)Tj     J    ^    J.f         J.o    J 


c     =  - 


For  the  present  approximate  calculations  we  shall  ignore  the 
production  of  bromine  and  hydrogen  atonns  at  the  equilibriu.ri  conditions,* 
In  this  case  it  is  found  from  Eqs,     (12)  and  (13)  that 


♦  This  assumption  is  in  line  with  the  steady-state  postulate  but  is  not 
correct  at  elevated  temperatures  where  appreciable  concentrations  of 
bromine  atoms  may  form  in  bromine-rich  nnixtureso    The  required 
corrections  have  been  worked  out  and  have  been  found  to  be  relatively 
small. 
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(Yj-Y^  X   ^  (Y2-Y2  f)h° 


(Xj-Xj  ^)Wjh°+(X2-X2  f)W2h° 


^    I,£      1,0'     1  1      ^    2,i      2,0'     2  2 


since  h^  s  0.    But 


and 


^^2  "  ^2,£^  ''^^r^i.f* 


<^2,f  "^2.o*=2<^l,f  "  ^1.0^ 


Hence 

P  * 


We  ahall  be  interested  only  in  flames  for  which  X,       -  0.     Hence 

1,0 

the  energy  equation  reduces  to  the  following  simple  expression: 

mc  1,  f 

P 

Following  the  conventional  procedure  for  solving  problems  in 
laminar  flame  propagation,  we  nov/  divide  the  continuity  equation 
given  in  Eq.  (U)  by  the  energy  equation  given  in  Eq.  (14)  in  order 
to  eliminate  the  distance  x  as  the  independent  variable.     In  this 
manner  it  is  foxmd  that 
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m^     dX,  k,    ^^2  3/2       (a--J)(l-a~ji) 


3/2 


l-a4.^)xJ[(9-l).(1.0^)(i.X,/X,^,i] 
We  now  introduce  the  assumption^  '  that 


.   (15) 


1/2 

1/2   -0,/0 
4  =  ^2®        « 


(16) 
(17) 


kj      Bj      -Oj/0 
*5     **5 


(18) 


Equation  (15)  now  becomes 


dXj 


A,W     B, 

m  c  5 


1/2 


^Z^^,^ 


3/2   -1/2     -©^/9 


(19) 


X    3/2 
(2a-X^)(l-a-^) 

X 


^l.ai4  .  ^)xj  [(0.1).(1.QJ(1.^  j 


wnere 


«'  =  e,  fie. 


(20) 


'2  -  ^2  ^  2  "1       • 

It  is  now  convenient  to  introduce  the  parameter  ^i    through  the 


relation 


m  c 
P 


1/2  3/2 

5  ^^f 


(21) 
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whence  Eq,  (19)  becomes 

X      ^^^ 


(22) 


The  boundary  conditions  are: 


d  -  0, 


when 


Xj  =  0 


(23) 


and 


e    :::    1 


when 


X.    -  s  2a        for  a  <  l/2 
X.    -  a  2a        for  a  =  l/2 
J   J  =:  2(1  -  a)  for  a  ?  i/2 


S-    (24) 


X 


The  solution  of  the  boundary  value  problems  is  discussed  in  the 
following  Sections  III,  IV,  and  V  for  a  <  1/2,  a   =1/2,  and  a  y  l/2 
respectively. 


/ 
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:•!.     BROMINE -RICH  FLAMES 

For  the  bromine -rich  flames,  solution  of  the  boundary -value 
problem  can  be  effected  conveniently  by  use  of  the  von  Karman- 
Millan  technique.      '    The  boundary  conditions  are 


and 


0  =  8j       when  X,  =  0     1 


0  =  1        when  X.  =  2a 


J 


From  Eq.  (22) 

^     -0^/0  ^  ^/2 


i^r- X 


|0-l)i(l-0^)(l-:^)         (l-a-0381X^) 


(24) 


(2S) 


where  the  appropriate  numerical  value  of  k./k^  has  been  introduced 
;i.3e  Section  VIC  for  details). 

Rearranging  Eq.  (25),  we  get 


.*.      I      (l-a-yX. ) 


1^1 


(©.1),<1-0^)(1.^) 


I^ 


2a -X 


1 


f roxn  which  it  follows  that 


71. 


-0' 
''2 


u      /d®    X      (1-1. 762a) 
lim  (h«-)  =  -S 4/2^ 

©-♦1  «^i         (l-2a)'^''^ 


-I 


(26) 


(27) 


or 

lim     ,  dO  , 


Ae''^    .iLlLi§^ 


(l-2a) 


572" 


.  ^^•^o^I-1.762a) 


TiT 


(l-2a) 


Td 


(28) 


Defining 


5  (a) 


1.1.762a 


(l-2a) 


TTT 


(29) 
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and 


1-9 


b'(a,0^)=-2: 


(30) 


then  Eq.  (28)  becomes 


e--*.i<35r>' ::^ 

*       Ac     ''  +<5{a) 


(31) 


Introduction  of  Eq.  (31)  on  the  left  hand  side  of  Zq.   (26)  for 
0^1,  X./v/  2a,  leads  to  the  following  relation: 


-0 


-4.C 


2  b'(a,Q.)        ,     (9-i).(l-9^)(l-^) 


Jle         >6(a) 


2a-X, 


or,  using  Eq.   (30), 

1-e  =  b'(a,0^)(2a-Xj) 


-0» 


1 -ZKT- 

-/le 


5(a) 


or 


1-9  =  b'(a.9^)(2a-Xj) 


6(a) 


— igT 

Ae     ^+*(a) 


(32) 


Introduction  of  Eq.  (32)  on  the  right-hand  side  of  Eq.  (26) 


yields 


A9 


i.     -9'  /O 

?       ^2^*'d9 


(l-a-O.  331Xj) 
(l-a-^X^)^^^ 


b'(a,9^) 


6(a) 


- — rg'r 

Ae      ^    5(a: 


;(33) 


integration    of  this  last  relation  between  9^,  X.  and  9  =  1,  X.  =  2a, 

2a 


leads  to  the  result 
I       1 


\    9"^e"'^^'d9=^l§J«a) 


-0* 


JLq        +i(a) 


(l-a-0.3fcIXj 
J ^  dX,  .     (34) 

(1-a-^  xy'^       ^ 
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For  9^9^,  Xj'vr  0,  we  get  from  Eq.  (26) 


.'1  '^y^  do 


(0-l)+(l-Q^) 


2a 


(l.a.0.331Xj) 


(35) 


the  integral  of  Eq.   (35)  from  0^ ,  0  to  0. ,  X    is 
®t  1 


0  "e 


«j  o 


do  = 


la  A. 


o 


(l-a--ir  X,) 


7^1 


But 


(8) 


0. 


dO  ^ 


(^TTFrr:^ 


r*  4  -®2/« 

J   e  ^  c  '^  de 


whence 


C^  4   -®2/® 

J   ©  ^  e  ^  de- 


2a -A 


*(l-a-0.381 
i        (l-a-ixj 


X,) 

^dX,.  (37) 


do 
Sinc« -jw-  is  to  be  continuous  at  9^*^*.  ^^  follows  from  Eqa. 

(35)  and  (33)  ttiat 


(O^-DHI-Qq) 
Za 


-1    -d'/6. 


(l-a-0.381X^) 


AO.   7e      ^     ^     Cl-^-iXj) 


whence 


l-a-0.381X 


'       n-W-  '^'<^'«o> 


t 


l- 


llal 


-y: 


(o^-D+d-e^) 


2a  ^ 


-0' 
b'(a.e^)e     ^ 

Tie        +5(a) 


Tie      "icS(a)^ 


(38) 
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In  view  of  Eq.  (38),  Eqs.  (34)  and  (37)  can  be  added  to  yield 
the  result 


'9*2      rZB. 


-^    -9«/e  b'(a.O)c     ^      \       (l-a-0.381X 

J      ,    ^e     2     dQ= :4 J  ^         .jl     dX^  .  (39) 

0,  ^e     ^  +  S(a)       0     Cl-a-jXj) 


Xj) 


In  order  to  evaluate  the  integral 
1 


\ 


4  -«2/® 

e    '^  e      "^     dO 


we  set 


I  121  4.,^  2, 
=  y.    u  =  y  .  — T-  =  y      and  dQ  =  -  -^  dy  . 

TT  ^  y 


Substitution  of  these  relations  yields 


J".-^ 


.©•y2 


I 


1 


-0^y 


7     -e^/e 

dd  =  -2         I       ye  -i-  dy  =  2  \        e      ^     dy. 


But 


)^  , 


2. 


-9Jy2 


ix. 


e  dy  =  J        -20'  ye 

I  1  -2V92 


or 


.e'y2 


2yo; 


15 


-*i 


2    J      y 
1 


Comparison  of  the  two  preceding  relations  leads  to  the  result 


J      S    ^  e     ^     dO  =  -40^^  \ 


-13- 

e  dy    + 


1 


-e- y2 


?55r^ 


JL 
1 


■i^ 


'     -^2^',         f>F        -«^2y' 


-«2,^.-^2/«i 


40*  e     *"      dy-\'*'i      e     *■      dy     +2e     ^-Z  W^ 


Q 


-m 


I     -0'  -Oi/Q.\ 

=  (2e    "^  -2  yEr  c   ^    ^)  ^^Y^ 


e"*  dt 


dt      . 


For  sufficiently  large  values^   '  of  ^U^i  ^®  obtain* 
1        1 


e    '^  c     "^      dO  =     2e     ^ 


2  \e"^  dt-2\  e°^  dt 


7  _    — 2'-i 


=    2e     "^  -  2  yS^e 


-e;/o. 


-2)^ 


1  -E2()n5T) 


] 


where 


<T>' 


But 


i-E2(>ej) 


-0' 
*'2 


J  zy^ 


1        1       ^    1.3         1.3.5    , 
1  -  ygr-   + , 5  !■  . . . 


*Xn  the  practical  applications  considered  here,  01 /O.  is  sufficiently 
large  to  justify  the  approximations  made  here. 
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whence,  finally. 


r 


-^   .9'/©  -©•        -e' 

Q    ^  e     ^      dOcfZe     ^  -2e     "^ 


1         ^      h    ^'^ 


l-3'5 


f  . 


or 
•1 


-2  '^y^ 

0         e 


•9 


dO^     ^ 


1  -    j^-3     ,  1-3' 5    ^ 
^^^2^       (20p^ 


(40) 


71, 


-a* 
*'2 


-a' 
2 

Equation  (39)  can  zu>w  be  solved  for  Jle         with  the  result 

a-0.381XildXi 
. — ^ ;__* 


(l.d^)e 


-82  ^2a 


2a 


1«3    ,    l-3;^5 

^2^ 


(^^      (20')^ 


l«1.762a 


or 


-03 


TLe 


(i-9,)e; 


J     (1-a-O. 


331  Xj) 


dX, 


0     {l-a-i/S!lti)3/i;^'^l 


(1-1. 762a) 
(l-2a)^^^ 


(4i) 


where   d(a)  and  b'(a,  0  )  have  been  replaced  according  to  Eqe.  (29) 
and  (30)  respectively.    It  is  now  necessary  to  evaluate  the  integral 


2a  /^2a  /-2a 

f       (l-a-0.381X.)  \  ...  \         X   dX^ 

0       (l-a--jX^)  Q     (l-a-^X^)  Q     11-a-^Aj^ 


jn 


Evidently 


-lis. 


r 


, ^1^  =4(i-a)    -—-—:: 


)m 


and 


2a 


5X,dX 


0  (i-^-7^i) 


S7I 


=  0.381 


8a 


i  i6|y  i-2a  -rr^ 


Therefore 


2a 


(l.a-0.38ixp  J 


4(l-a)-3.  05a-6. 10fl2. 20a 
-i— r  4(l-a)-6. 10(l-a)l 

m^  I  J 


^  5.15a>2.1Q    ,  2.J^l-a)^  ^^^j^  ^^^^ 


The  fiaal  result  for  bromine -rich  flames  becomes 


Ae 


-oj       (i.e^)eif(a) 


2a 


1- 


1-3   .  1-3-5 

2^ 


^^  (2en 


(U1.762a) 

*     3/^ 

(l-2a)''''^ 


(43) 


froin  which  the  eigenvalue  A  can  be  obtained  for  bromine -rich  misi- 
tures.     The  mass  burning  rate  m  is  then  determined  through  T^<c\,  {Zl). 


-lo- 
IV.     THE  STOICHIOMETRIC  MIXTURE 

For  the  caae  a  s  v  and  X.   .  =  i  the  banlc  differential  equation, 
£q.  (22),  becomes 


1      .guO  5/2 


or 


^  '"^-"'^  *    (I'-'a^viik^)   ^  »(9-i»(i-v«»-V  •       «^> 


In  order  to  integrate  £q.  (44)  we  apply  an  ingenious*  as  yet  un> 
published,  method  invented  by  Dr.   Theodore  von  Karman.      We  replace 
the  singular  point  (0  »  1,  X^  =  1)  by  (0  «  1,  X^  a  XJ<  1);  then  Zq,  (44) 
becomes,  for  (0  -  1)  small  compared  to  (1-©  )(l-X.), 

Integrating  £q.  (45)  from  (O^.Xjt*)  to  (l,xp  leads  to  the  result 

0*  x^* 

For  6*/^  1,  the  preceding  expression  becomes 

XJ 

Jl      "®2  C         (1-0.761  X.) 

-re        (l.©^)o:;(l-e  )      \       -_1    dX     .  (46a) 

US  °     J        (i-xy'"^  ^ 

Xf 

If  Eq.  (45)  holds,  as  a  first  approsdiriation,  from  (9^,0)  to 
(1,X|).  then 


-17- 


Ji      e 


-9;/e 


/7j  y^ 

0. 


do  =(1-0^) 


^1 

f      (1-0.761X^) 


dX, 


(47) 


or,  utilizing  the  result  obtained  previously  for  the  integral  appearing 
on  the  left  side  of  Eq.  (47), 


-Q! 


vl      e 


1^  IT" 


1- 


1-3        1-3-5 


^^       (26^) 


■^^^     •     •    •    •    • 


XJ 

[  (i-a7 

°    ^        (1- 


-a761X,)dX 


0       <'-''i> 


From  the  original  differential  equation  it  is  apparent  that  the 
limiting  slope  at  the  hot  boundary  is 


^  =(i-e  ) 


to  which  we  join  the  curve  of  £q.  (46a)  at  (9'^,X^).     Therefore,  at 
(0*.Xf). 


(1-0*)  =  (1-0^)(1-X*)   . 


(19) 


From  I  qs.  (46a)  and  (49) 

-0 

-     € 

dO 


^   e   ^2  ^j 


^1 

f       (1-0.761X,) 

-^1*  >  =   \    rrr  ^^i 


^1* 


WithgJ  =  (l-®o).   Eq-  (45)  becomes  at  (0*,  X*^): 


vt 


-01 


3/2 


^     (1-0.761^^^*)    <^"^t^ 


=  1. 


(5X)) 


(51) 


We  next  integrate  £q.  (50)  noting  that  Xf^  and  X|  are  very 


close  together  so  that  we  may  approximate  (1  -  0.761X.)  in  the  inter- 
val Xf  to  XJ  by  (I  -  0.761X'J).     Hence  Eq.  (50)  becomes 


A      -»J 


1^   e      "  (l-X*)  =  (1-0.761X 


C         dX 

f)    \    ^7p  =  2(1-0.  761X*) 

*     cJ     (l-X,)^'*^  * 

X, 


/(1-X')(1-X«J) 


or 


-0'  3/2 


/rxfyrixj- 


=  2 


yiTsn 


(52) 


From  Eqs.  (51)  amd  (32)  we  now  see  that 


or 


-2  =  1 


/ 


TTxY 


3 
7 


whence 


xi4xj-| 


(53) 


The  next  step  is  to  obtain  an  explicit  esqsression  forXj  frora 
Eqs.  (48)  and  (50).     Equation  (50)  nnay  be  written  in  the  form 


-«2 


XJ 


^       e .  1  C  (1-0.761  Xi)dXi 

X«| 


(54) 


or,  using  Eq.  (48),    ,j., 

(1-0   )       r    Ml-0.761Xj)dX^ 


37T 


1  - 


1-3 


1-3-5 


i^)^     (2^ 


(l-C.761X^)dXj 

jn — 


9^(l-Xf)     J  (l-Xj) 

X>l| 


_f  n_ 


Replacing  X*  in  terms  of  Xi  according  to  Eq.   (53)  in  the  preceding 
expression  leads  to  the  following  explicit  relation  for  Xh 

[1  ^1  "^1 

^'zO^^iZOn^'  '"]  \  (1-0.761X,)  \(1.0.76U9 

1^1 -I  ' 


Integration  of  Eq.  (55)  leads  to  the  following  first  approximation 

forXJ: 

fl     1-3        l'3«5  ] 

xj  «  1 — — •: — ' .  (56) 

Hencet  from  £q.  (53), 

^t-^ 3d'  (1-6 J .  (57) 


and, from  Eq.  (49)» 

(,1-3      ,    1-3-5  "I 

,,^,,    i''m)'mj^"-i  (58) 

>  30' 

Comparison  of  Eqs.  (56)  and  (57)  shows  that  X|  and  X^  are  indeed 
so  close  together  as  to  justify  the  approximation  made  in  the  deriva- 
tion of  Eq.  (52). 

-0' 
2 
In  order  to  obtain  a  first  approximation  to  7l  e         we  proceed 

from  Eq.   (54)  as  follows: 

^      -0^       2(1-0. 76lXj*)       [        J 

,9 


l-0.761X^  1-0.761(|-X|  -|) 


after  replacing  X|  according  to  Eq.  (56)  it  is  found  that 


n 


0.239+0.761 


,     1-3     ^  1»3'5  . 
(2Q^)      (2e!>)^ 


5(6'   (1-6   J 


fl     1*3     .    1-3' 5  1 


T72    • 


(59) 


38^  (1- V 

Therefore*  for  the  range  of  values  of  Oi  actually  encoiintered, 
Tte     2 


C:"    1"^     >l^\{\'^^    [  0.761  +  0.717  e2<^-®o^  1 

from  which  the  eigenvalue  -A,  for  the  atoichiometric  mixture  ratio 
can  be  calculated. 


(59  a) 


V.     HYDROGEN-RICH  MIXTURES 

The  mathematical  technique  which  we  have  employed  for 
stoichiometric  mixtures  can  be  used  also  to  solve  the  boundary- 
value  problem  for  hydrogen- rich  mixtures.     However,   since  the 
calculation  is  somewhat  laborious,  we  follow  another  suggestion 
by  Dr.  von  K^rm£n,  according  to  which  we  test  a  more  highly  ap- 
proximate procedure  for  nriixtures  containing  43.  3°/o  Br^*  for 
which  the  solutions  become  particularly  simple.     The  success  of 
the  simplified  procedure  for  mixtures  containing  43.  3  /o  Br^  gives 
us  confidence  in  the  results  obtained  for  hydrogen- rich  mixtures, 
in  general,  for  which  we  utilize  a  method  of  calculation  in  which 
the  limiting  slope  of  the  0  vs.  X.  curve  does  not  have  the  correct 
value  at  the  hot  boundary. 

A.     The  Mixture  Containing  43.  3  /o  of  Bromine 

For  mixtures  containing  43.  3  /o  of  Br^  a  considerable 
simplification  is  noted  in  the  basic  differential  equation.  For 
tliese  mixtures  X.  ^  s  2(l-a)  and  it  is  seen  that 

*      2a  =  ^'^^y     .       l-a  =  0.433; 
Equation  (22}  becomes,  for  (0-1)  small  compared  to  (1-0   H^'^i)* 

A  iflt —       ®  ®  55r  ~  •  ^60) 

0.3K1  ax^  ^,^j 

(1)    Direct  Integration  of  Eg.  (60) 

Without  regard  to  the  li]::[iiting  slope  at  the  hot  boxindary, 
we  can  integrate  Eq.  (60)  directly  with  the  result 


VI 


2(l-a) 


J — m*^® srnr  J 


«i 


dX, 


STTTI 

(l-a-ji) 


(61) 


since  1-a  3  0.433.    Hence 


Tie 


-0' 


•2 


and 


1 .  i'3  )  r3»5 


.  2(l.a) 


0.381(l.d^)(- 


inrsT 


m  2.318(1-0^) 


^/: 


M.e 


-0^        2.318(1-0^)0^ 


^     1-3     ^   1«3'5 
^^®2^        (20')^ 


(62) 


The  fact  that>^£q.  (62)  yields  a  fair  approxin^tlon  to  the  solution 
of  the  boiindsiiy-value  problem  will  now  be  established  by  a  more 
careful  integration  in  which  dO/dX.  has  its  correct  limiting  value 
at  the  hot  boundary.    On  the  basis  of  this  result,  we  shall  inte- 
grate the  flame  equation  for  all  hydrogen- rich  mixtures  directly 
in  Section  VB. 

(2)    Integration  of  Eg.  (60)  With  the  Correct  I-.imiting 
Slope  at  the  Hot  Boundary 

Integration  of  Eq.  (60)  from  {9*,  XI)  to  (1,  Xl)  leads  to 
the  relation: 
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.9'/0 


,0.381, 


^  \  ~^<^«  =  <^-V<^^) 


dX 


1 


(1-a-.^) 


^frr/2 


0* 


Xf 


(63) 


For  0*  ^  1 


-0' 


^1 


-Ac     **  (1-0*)  =  0.880(1-0^) 


dX, 


(l-a-  ^ 


!S7r7Z 


A 


If  £q.   (63}  holds  from  (9..  0)  to  (1.  xp»  then 
I  X' 

dX 

de  =  0.880( 


i 

J        172" 


X' 


d^     0 


■xp72 


0    (l-a-  -ji  ) 


(64) 


or 


A 


-0* 

7 


^2 


*1 


1.  ^'^  +illll- 


dX, 


^  0.880(1-0^) 


(l-a.-5~) 


xrrrz  • 


We  now  join  S,<i*  (64)  to 


^  1-0 

do  ^  ''o 


S^^  "2ITiiT 


at     (0^,  X-J) 


(65) 


since  dO/dX^  »  (UO  )/2(l-a)  is  the  correct  limiting  slope  obtained 
from  Eq.   (22).    At  0  =  0*.   X,   =  X'J: 


1-0 
1-0*  *      o 

2(l.a).Ji'{    'I(T:^ 


(66) 


and,  from  Eqs.  (64)  and  (65), 
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TLe 


(1-a) 


dX, 


0. 880  \ 


3rTr/2 


With 


£q.  (60)  becomes  at  (©*,  X^)  for  0* 

x«^/2 

(l-a-y) 


^c 


"2 


2(l.a) 


s    0. 3S0 


(67) 


(68) 


From  Eq.  (67) 
-0 


Ai 


^^ff!iF^--^•"''  u->-5i''M<-4) 


.♦,  1/2 


.  (69) 


But,  from  Eq».  (68)  and  (69): 

.-  1/2 

)^  X  0.  880  [2(l*a)-Xf  J         =  3.  520 


r  ♦ 

xrTi/2 


|(l.a)-f  J         .[(l.a)-.ji 


XM/2 


or 


[(1-a). 


Xf 


il/2 


.il/2 


(l-a). 


i]         4[(i.a).^] 


♦,1/2 


whence 


x: 


(l.a-U)=|(l.a.9L)    . 


(70) 


From  Eq.  (67): 
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-0'  P 

-Ae     ^       2(I-a)0. 880  I 


Xf 


dXj 
i(l-a).J.l 


and,  using  Eq.  (65), 


0.880  (1-6   ) 


1  ,  1.3   ^1.3- 
^^®2^     (20! 


X' 


•y  r  ^^^i  0. 880(l-a)     f      *^^1 

TTTT ^  \   r ^TF^" — f — ^  \  f ^qm 

^  0  X^ 


Thus 
(1-a) 


1    .J'?.  1   ^'•^'? 


(l-»o>«2 


F^j^TT 


^f 


1/2 


vM  1/2 


{l-a-  -w-  )       -  (1-a-  -*-) 


1/2  X 


,1/2 


(1-a)        -  (1-a-  -»-) 


Using  Eq.  (70)  to  eliminate  X|,  the  preceding  expression  becomes 


I  fi      1*3        1*3*5 

|(l-a)[^-(:SBJ)^(15Ty?-- 

(1- V  «2  p-^>  -  -r  J 


(1-a) 


1/2      1 


X 


«  1/2 


(1-a-  -^  ) 


which  yields 

Xt  ^^^  1/2 

(l-a-^M        =(l-a)        <1-. 


1  - 


tt^5;t^2 


1-  ^'^    ,  -  -  : 


l-3-5|         I 

I 


•   (71) 


Combining  £qs.  (68)  and  (71),  and  noting  that  (I -a)  »  0.433,  we  get 


Ai 


-0' 


1.16 


I   - 


1 


v^^'^r^z 


.      1'3         1»3'5 


(29')^ 


I 


/2 


(72) 
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The  fact  tkat  Eqs.  (62)  and  (72)  yield  simila  '  results  can 
be  seen  by  expanding  the  denoxninator  of  Eq.  (72)  for  large  values 
of  (1-0  )  0^.     In  this  manner  it  is  found  that 

Jle     ^  :i   2.32(l-©^)9^  (72a) 

which  is  seen  to  be  practicadly  identical  with  Eq.  (62)  for  large 
values  of  OL.    On  the  basis  of  this  result  we  are  encouraged  to 
integrate  the  basic  differential  equation  for  hydrogen-rich  mix- 
tures directly  between  0. ,  0  and  1,2(1 -a)  after  neglecting  the 
ternn  in  (0-1)*     Details  concerning  this  work  are  described  in 
the  following  Section  VB. 

V«B.     Hydrogen-Rich  Flames  with  Arbitrary  Concent  'ations  of 
Bromine 

The  ^result  for  the  general  case  of  hydrogen-rich  flames 
is  obtained  most  conveniently  by  using  direct  integration  of  the 
differential  equation  with  the  termi  in  (d-1)  neglected.     For 
hydrogen  -  rich  flaznes    a>T  and  X,  ^  s  2(l.a).     Equation  (23) 
becomes 


J  gl/2  •  »-*    %  (2a.X^)  [2(l-a).xJ 


But 


-■'7™ 


2(l-a)  2a 

o      7r~rTT:r~mT7r  =  J    — 7=; 


762a 


dy 


0  (2a.xp  [2(i.a).xJ^      ^^.^ 


yyy-(4a-2) 
2a 


or 


2tl-a) 
(2a-xp  [2(l.a).X^j^^^ 


[1-1.7^ 

=  [ — ^rn 


-I. 762a 


Z:i 


COS 


-ll/JaTT 


2y7^N2a 


and  £q.  (73}  becomes 


-1   -0^/e 


do  s 


)^ 


(1-0  J 
fl   X  0.381     ,^J\ 


1-1. 762a       2       ,„-l)grT  .  ,v7-7r        /7a\ 

6.S81   ^^''''^    [-^j+2/irzi|.(74) 


Since  0    is  not  a  sensitive  ftinction  of  a  we  shall  treat  it  as 


a  constant.     Then 


^     -O'/O 


^  J    -^i dO  =  0.  539(1-0^)  g(a) 

0      y^  ^ 


(74a) 


where  g(a)  is  defined  as 


g(a) 


=    1 


rri* 


l>1.762a 

(5.3SI 


/iim 


cos 


t/^ 


+  2  KITSI 


(75) 


Finally,  integration  of  Eq.  (74a)  yields  the  result 


JL  e'®^  c:      ''  i^V'" Va!^  ^^^^   .  (76) 


f-  (1%)   ' 


2'      (20^ 

from  which  the  eigenvalue  -A.  for  hydrogen>rich  mixtures  can  be 

calculated.     For  (1-a)  =  0.433,   1.762a  =  1,  and  Eqs.  (75)  and  (76) 

reduce,  of  course,  to  Eq.  (62). 

-0' 
We  have  now  obtained  analytic  expressions  for  Tie         for 

all  possible  hydrogen  ^bromine  mixtures.    Before  summarizing  the 

basic  expressions  for  the  laminar  burning  velocity,  we  proceed  to 

consider  the  relevant  physico-cheznical  parameters  which  must 

be  used  in  order  to  obtain  numerical  results. 
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VI.     PHYSICO-CHEMICAL  PARAMETERS 

A.     Equilibrium  Constant  for  the  Reaction  Br,^!^  2Br 

The  value  of  the  equilibriiurn  constant  K    (i.e.  ,  of  the  equilib- 
rium constant  expressed  in  terms  of  partial  pressure  ratios)  for  the 

bromine  molecule -bromine  atom  system  have  been  taken  from  the 

(13) 
Bureau  of  Standards  compilation^      '  which  covers  the  temperature 

range  between  298^K  and  ISOO^^K.     Equilibrium  constants  at  1600^K 
and  1700  K  were  obtained  by  extrapolation  of  the  Bureau  of  Standards 
data.     Equilibrium  constants  as  a  function  of  temperature  are  sum- 
marized in  Table  I. 

The  data  given  in  Table  I  have  been  fitted  to  an  expression 
of  the  form 

by  plotting  log(K   /T)  vs.  (1 /T).     The  plot  is  s/iown  in  Fig.   1.    Ap- 

It 

propriate  nuxnerical  values  derived  from  Fig.   1  are 

A  s  45,  210  cal.  /mole        and  B  »  230.  5  atmos.  /^K  . 

The  ejjuilibrium  constant  K  »  expressed  in  terms  of  concen- 

••*  c 

tration  ratios,  for  the  reaction 

Br2=5=^2Br, 
is  related  to  K    through  the  expression        K    »  K   /RT. 

Hence  it  follows  that 

V        230. 5T         -45,210/RT 
*^c  °   82.07T     ^ 

-i   ono     -45,210/RT  ,       ,       3 

=  2.308  6  moles /cm 

if  concentrations  are  expressed  in  moles /cm  .    But 
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where  k.  and  k^  are  the  specific  reaction  rate  constants  for  the  pro* 
cesses  described  in  Eqs.  (1)  and  (5),  respectively. 

In  terms  of  the  parameters  needed  for  the  calculation  of  lami- 
nar burning  velocities  we  find,  therefore,  that 

Bj/Bg  a  2.808  moles /cm  (77) 

and 

0^    =45,2lO/RTj.  (78) 

B.     The  Specific  Reaction-Rate  Constant  for  the  Reaction  Between 
Br  and  H. 


'2   ■ 


The  specific  reaction-rate  constant  for  the  process 

^2 

Br  +H2— =-^   HBr  +  H 

is  based  on  careful  e9q>erimental  studies  carried  out  by  Bodenstein 

{14) 

and  L,ind.  ^     '    The  experimental  results  have  been  re -evaluated  by 

(12) 
Pease^     '  who  finds  that 

*k2  =  4.  56  X  10^  Vr'e-^®'"^®^/^'^  cm^ /(mole-sec), 
•-  (79) 

whence 

B,  =  4.  56  X  10^*^^  yrr       cm^/(molc.sec). 

We  shall  utilize  Eq.  (79)  for  the  specific  reaction-rate  constant  k. 
in  preference  to  slightly  different  forms  i;  roposed  by  other  investi- 
gators.^  *    '    However,  it  is  important  to  note  that  alternate  correla- 
tions of  the  available  experinnental  data  lead  to  greatly  different  pre- 
dicted values  for  k,  at  elevated  temperatures.     Thus  Hirschfelder 
and  his  collaborators'    '  recommend  the  relation 

k,  =  3.45  X  10^°  T  c-^^'^^^^^'^cm^ /(mole-sec).  (79a) 
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A  comparison  of  the  calculated  values  of  k^  determined  from  Eqs. 
(79)  and  (79a)  Is  aho>vn  in  Table  II  for  the  temperature  range  which 
is  of  interest  in  connection  with  calculations  of  laminar  burning  ve- 
locities for  H^-Br,  mixtures.     Reference  to  Table  II  shows,  for 
example,  that  the  values  of  k2  calcxilated  from  Eq.  (79)  at  ISOO^K 
could  be  too  large  by  roughly  a  factor  of  1.  7.     Hence  the  laminar 
burning  velocities,  which  vary  as  the  square  root  of  k,,  could  be 
too  large  by  perhaps  a  factor  of  1.  3. 

C.     Ratio  of  the  Reaction  Rates  ^4/1^3  for  the  Removal  of  Hy- 
drogen  Atoms 

The  ratio  ox  the  reaction  rates  k. /k>  for  the  removal  of  H- 

atonas  according  to  the  processes 

^^3 
H+Br2  — ^  HBr  +  Br 

and  k. 

H+HBr—^H^  ^Br 

is  again  based  on  the  experimental  work  of  Bodenstein  and  Lind^    " 


and  Bodenstein  and  Jung.  ^      '    The  value 


k^/kj:»  1/8.4  (80) 

does  not  vary  with  temperature  over  the  relatively  narrow  temper- 
ature range  for  which  experimental  measurements  are  made.     We 
assume  that  this  result  can  be  extrapolated  fro'in  about  600  K  to 
temperatures  as  high  as  1700^K.    Since  the  activation  energies  for 
the  reactions  (3)  and  (4)  are  probably  very  small,  this  extrapolation 
should  not  introduce  a  large  error  even  if  the  activation  energies 
for  reactions  (3)  and  (4)  are  not  identical. 
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D.  The  Adiabatic  Flame  Temperatures  and  Equilibriuin  Composi- 
tions for  H2-Br2  Mixtures 

Adiabatic  flame  temperatures  and  equilibrium  compositioas 
have  been  determined  by  the  use  of  standard  procedures.       '    The 
results  for  T,  as  a  function  of  composition  of  the  gas  mixture  are 
plotted  in  Fig.   2. ;  the  equilibrivun  mole  fractions  of  HTr,  H,,  Br,, 
and  5r  as  a  function  of  composition  are  shown  in  Fig.   3.     Relevant 
numerical  values  are  summarized  in  Table  III. 

E.  Thermal  Conductivities  in  H^-Br^  Misstures 

Values  of  ^^  as  a  function  of  composition  for  H^-Br^  mix- 

'  (5) 

tures  have  been  calculated  by  Cooley  and  Anderson^   '  on  the  basis 

of  the  elementary  kinetic  theory  of  gases,  which  does  not  include 
proper  allowance  for  the  influence  on  thermal  conductivity  of  energy 
transfer  between  internal  and  translational  degrees  of  freedom.     The 
results  of  Cooley  and  Anderson,  which  represent  lower  limits  for 
the  thermal  conductivities,  are  summarized  in  Fig.  4,  where  we 
have  plotted  "A.£  as  a  function  of  composition.     Reference  to  Fig.  4 
shows  that  ^>increase8  slightly  with  T£  and  more  rapidly  with  the 
concentration  of  H2  in  the  gas  miaeture. 

Actually  the  problem  of  estimating  appropriate  therjcnal  con- 
duction coefficients  is  greatly  complicated  for  flames  by  the  inevi- 
table con^position  clmnges  occurring  during  reaction.     Campbell  and 
Hirschfelder^   '  propose  a  relation  of  the  form 

A»A'-s-?l"T.  (81) 

where  both  ^'  and   (d"  are  parameters  which  must  be  determ.ined 
separately  for  each  gas  mixture.    Our  formulation  of  the  problem 
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of  flaiTse  propagation  does  not  permit  the  use  of  this  two -parameter 
expression  for  ^  without  a  complete  recalculation  of  the  results, 
because  of  the  change  in  the  assumed  temperature  dependence  of  A. 
In  view  of  the  crudeness  of  the  approximations  for  A  ,  we  do  not 

feel  that  the  additional  labor  is  warranted.  For  this  reason  we  have 

lid) 

used  the  following  relation         for  the  thermal  conductivity  of  a 


mixture  of  species  i  and  j: 


•^f  =='^mixture<'^f>=I 


X  .X,+  A  ,X.  f ~ — ^ , 


<82> 


Equation  (82)  yields  reasonable  results  for  binary  mixtures  in  which 
the  pure  components  have  greatly  different  thermal  conductivities. 
Thus  £q.  (82)  is  particularly  v/ell  suited  for  making  estimates  of  A  £ 
in  hydrogen-rich  flames.    In  order  to  calculate  the  thermal  conduc- 
tivities of  the  pure  components,  we  utilize  the  following  relations 

(6) 
proposed  by  Can^pbell  and  HirscMelder*     '  and  neglect  again  the  con- 


tributions of^Br  and  H  at  T£  to  the  thermal  conductivity: 

^    AH2(T£)  =  5. 13  X  10"^    /-5^    .  (83) 


,     -P^ — 
^Br2(T£)  =  1.20x10-^^^     ,  (84) 

6.  OS  X  10'^  fTZ 
^HBr(T£)  = i^ (85) 

f  o 
where  Af  ^b  expressed  in  cal/(cm  sec     K). 

Numerical  values  for  A  r  as  a  function  of  mixture  ratio  are  given  in 

Table  IV  and  are  plotted  in  Fig.  4  as  a  function  of  composition. 

It  should  be  noted  that  our  method  for  calcxxlating  /\r 
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emphasizes  the  nuir.erical  values  at  temperatures  close  to  T,.     This 
procedure  receives  partial  support  through  the  well-known  result  that 
most  of  the  chexnical  reactions  in  a  flame  occur  at  ten:iperature6 
close  to  T,. 

F.     The  Average  Specific  Heat 

A  more  exact  expression  for  the  average  specific  heat  than  was 
used  in  the  derivation  of  £q.   (14)  can  be  obtained  from  Eq.   (13),  viz. , 


(85) 
W4  o      W5 

+  — ^X.    .h°  +— i   X. 


.i-l] 


since  X,        ~  X.       =  X^       ~  0  and  h,  =  0  if  the  subscripts  4  and  5  iden- 
1,0         4»o         5,0  3  *^ 

tify,   respectively,   the  species  H  and  Br.     Although  the  contributions 
of  H  and  Br  are  small,  they  have  been  included  in  Eq.   (85).     Average 
specific  heats  as  a  function  of  composition  have  been  calculated  both 
with  and  without  allowance  for  dissociation  into  hydrogen  and  bronnine 
atoms.     The  results  are  listed  in  Table  V.     Reference  to  the  data  given 
in  Table  V  shows  that  the  influence  of  the  atomic  species  on  the  aver- 
age specific  heat  is  negligibly  sznall  except  for  near -stoichiometric 
miixtures. 

If  dissociation  into  atoms  is  neglected,  then  simple  explicit  re- 
lations can  be  obtained  for  the  average  specific  heats.     After  introduc- 
tion of  appropriate  numterical  values'^  the  following  relations  are  ob- 
tained: 

for  X.   £  =  2a,   X,  £  -  X-       =  (l-2a)-(l-a)  =  -a,   i.  e. ,  for  a^  -y  , 

4 

—    _  2.  466  X  10  a .  .,.,. 

*=p- (15^.  3-157.  8a)Tf(l-6^)    '  <®^' 


**■  Appropriate  nunnerical  values  are 


Wj  hj=  -8660  cal/mole  and  ^2  h®   =  7340  cal/molc^^^^ 
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for  X,  £  =  2{i-a),  X,  f  *  ^^9  o  "  "i^*^)'  i.e.,  for  a  :^  -p 


T     -         2.466  X  10^  (l.a) 

since 

W=:  2.016a  •  159.8  (1-a)  =  159.8  -  157.8a. 

The  values  of  c"    calculated  from  Eqs.  (86)  and  (37)  are  plotted  In 
FI2.   5. 
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VII.     COMPARISON  OF  CALCULATED  AND  OBSERVED  BURNING 
VELOCITIES 

In  order  to  calculate  the  laminar  burning  velocity  as  a  func- 
tion of  mixture  ratio,   it  is  convenient  to  sumn^arize  first  the  rele- 
vant theoretical  relations. 

The  laminar  burning  velocity  u^^is  given  in  terms  of  the  param^ 
etex  yL  through  the  basic  relation  expressed  in  Eq.   (21),  viz. , 

p 

Introduction  of  appropriate  numerical  values  for  E  y,  B. /Bs»  &nd  R 
for  p  =  1  atmos.   leads  to  the  relation 

T 


./SI 


o 


u  =  S.  31  X  lO""  /—i — ^^—    Q  cm/sec  (89) 

p 

if  Wis  expressed  in  g/mole,  A  .  in  cal/cm  sec     K»  T     in    K,  Ic     in 

ir 

cal/g  °K.  and  T,  in  **K. 

For  the  eigenvalue  A  e         we  have  obtained  in  Sections  III, 

..9 

IV,  and  V  the  relations  listed  below. 
Bromine -Rich"  Mixtures  (a<  t»  ^i  f  "  ^a): 


^2  ^  ^  [  yr:^ fT^     i 


TLe        ~  -^ —^ — 1.3     .  1.3 


(43) 


(1-1. 762a) 

(1.2a)^^^ 

1 

Stoichiometric  Mixtures  (a  =  -y,   X,    ,  =  1): 

r^-^ — T-T-      1'3'5  -1^ 

.239  V  0.761    LI-    T7P^)^(^??rP   ""    ^ 

.   -e'       n  [  39' (1-ej 

Ae      2    2£  '-^-i ^ ^ 


36^71^ 


3/2 


(59) 
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-^   y  6  6^(1 -6^)  [  0.761  +  0.717d^(l-Q^)l     .  {59a) 

Hydrogen-Rich  Mixtures  for  2a  =  1 -a/0.  381,   X,   ^  =  2(1 -a): 

I    f^2>  {T^  ■•■J 

or,  in  somewhat  better  approximation, 
-0' 

^«    ^    = — ? 1 r.'^  1-3     .  l»7-5 TYm    ^"^^^ 

(2©')- 


Hydrogen-Rich  Mixtures  (a>y,   X,    ,  =  2-2a): 


®2<^-''o> 

-0'        0. 539     n      I  r     1    I    ^iL,  ,    n — r 

/I        2  (T-a)  I      i-1.762a  -L  /2a-l  »  ...,» 


%  +2  y2-2a 


* 


1- 


Explicit  relations  for  u^can  be  obtained  readily  by  combining 
Eq.   (89)  with  Sqs.   (43),  (59).  (72)  and  (76). 

JLaininar  burning  velocities  as  a  function  of  mixture  compo- 
sition, computed  by  use  of  the  preceding  relatione,  are  plotted  in 
Fig.   6  and  recorded  ih  Table  VI.     The  ;.ata  shown  in  Fig.   6  and  Table 
VI  correspond  to  the  two  extremie  sets  of  values  of  thertrial  conduc- 
tivity discussed  in  Section  VI.  £.     Also  plotted  in  Fig.   6  are  the  ex- 
perimentally deterzTtined  values  for  the  laminar  earning  velocity 
reported  by  Anderson  and  his  collaborators. 
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Reference  to  Fig.   6  shows  that  the  calculated  values  for  the 
burning  velocity  are  appreciably  larger  than  the  observed  results. 
Although  laminar  burning  velocities  determined  from  Bunsen  burner 
experiments  are  notably  unreliable,  it  appears  likely  that  the  major 
source  of  the  discrepamcy  between  theory  and  experiment  is  the  re- 
sult of  uncertainties  in  the  chemical  kinetics  data  and  in  the  thermal 
conductivities  used  in  our  calculations.    For  example,  on  the  basis 
of  the  ratios  of  rate  constants  listed  in  Table  II  •  it  is  not  unreason- 
able to  reduce  the  largest  calculated  values  of  the  burning  velocities 
by  as  much  as  a  factor  of  1.3,  thereby  producing  much  better  agree- 
ment between  theory  and  experiment  than  is  indicated  in  Fig.  6. 
Furthermore,  our  calculations  are  based  on  the  steady-state  approxi- 
nnation  for  hydrogen  and  bromine  atoms.     We  neglect,  therefore, 

the  influence  of  diffusion.     It  is  well-known  that  the  effect  of  diffu- 

(9) 
sion  on  burning  velocity  is  to  decrease  the  burning  velocity.^   ' 

Finally,  it  has  been  demonstrated  that  the  steady- state  approxima- 
tion for  bromine  atoms  cannot  apply  exactly  at  elevated  temperatures 
for  the  bromine -rich  and  near -stoichiometric  mixture  compositions. 
Details  concerning  the  corrections  required  because  of  failure  of 
the  steady-state  approximation  for  Br-  atoms  have  been  worked  out 
by  von  Karman  and  Millan  and  will  not  be  given  here. 

In  spite  of  the  discrepancies  noted  for  t'  e  absolute  values  of 
calculated  and  observed  burning  velocities,  it  is  evident  that  the  cal- 
culated burning  velocities  have  the  correct  functional  dependence  on 
composition.     This  result  gives  us  confidence  that,  except  for  the 
questionable  approximation  concerning  the  steady-state  approxima- 
tion for  Br-  atonns,  we  have  a  valid  first  approximation  to  the  physico- 
chemical  processes  occurring  in  a  hydrogen-bromine  flame. 
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Table  I,    I^ogarithm  of  the  equilibritim  constants  K  ,  for  the 
reaction  Br,  g±  2Br,  as  a  function  of  temper  - 
ature,  taken  from  reference  13 


T(  K)  -log  K    (p  in  atmos. ) 


298.16 
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28. 

318 

28. 

110 
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11. 

2754 
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8570 

7. 

0392 

5. 

6220 

4. 

4858 

3. 

5536 

2. 

7746 

2. 

1136 

1. 

5454 

1. 

0512 

■ 

5730 

• 

0175 

.42* 


Eqs.   (79) 

and 

(79a 

)  ae  a  function  of  temperature 

T<°K) 

■ 

k2  from  Eq.  (79)/k2  from  Eq.   (79a) 

500 

.69 

600 

.91 

1000 

1.42 

1500 

1.67 

1700 

1.72 
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Table  HI.     Adiabatic  flame  temperatures  and  equilibrium 
mole  fractions  aa  a  function  of  composition  (p  =  1  atmos. 
Tq  s  323     K,  a  3  initial  mole  fraction  of  H2) 


ill  ^^^   ^^  ^HBr,f  ^H^  f        ^Bt^  f  ^,f        ^Br,f 

.1  675  .20  .80  

.2  1014  .40  .60         8.0x10"^       

.3  1324  .60  .40  - 

.4  1585  .79               .21          1.6xl0"^  2. 5xl0"^    7   3x10"^ 

.5  1683  .95  2.7x10"^  2.9x10*^  2.2xl0"^    4.9x10"^ 

.57  1475  .86  3.5x10*'^     .085  3.3x10"^    7.9x10"^ 

.6  1395  .80  7.6x10'^     .162  4.3x10'^    6.0x10*^ 

.7  1182  .60  7.1x10"^     .378  1.2x10""^    3.4x10'^ 

.8  923  .40  .60  

.9  725  .20  .80  


—  .44- 

Table  IV.     Thermal  conductivities  Af  as  a  function  of 

conripoaition,  calculated  from  the  empirical  relation  given  in 

(6) 

Eg.  (82)  and  based  on  the  data  of  Cairipbell  and  Hirschfelder  ^   ' 

1-a  A^  xlO       (cal/cm  sec  °K) 


0.1 

43.2 

0.2 

39.6 

0.3 

33.1 

0.4 

25.0 

0.5 

14.1 

0.6 

8.8 

0.7 

6.0 

0.8 

3.6 

0.9 

2.3 
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Table  V.  Average  specific  heat  c*   afl  a  function  of  composition 

for  the  H^'Br^  flames 

l-a  c_(cal/g  ^K)(dissociation  neglected)*      c"  (cai/g  °K)** 

0. 1  0. 394  . 394 

0.2  0.213  .213 

0.3  0.151  .151 

0.4  0.119  .116 

0.433  0.116  .114 

0.5  0.112  .098 

0.6  '                      0.095  .076 

0.7  0.071  .065 

0. 8  0. 064  .  064 

0. 9  0. 043  .  043 


Recalculated  from  Eqs.  (86)  and  (37) 
♦  ♦Calculated  from  Eq.   (35) 
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Table  VI.     Calculated  laminar  burning  velocities,  vio^for 
Hg-Er,  flarnes  as  a  function  of  mixture  composition 


1-a 

0.1 

0.2 

0.3 

0.4 

0.433 

0.5 

0.6 

0.7 

0.8 

0.9 


f  55 
i^Based  on  the  thermal  conductivity  data  of  Cooley  and  Anderson.  ^   ' 

^'^Based  09  thermal  conductivities  calculated  from  Eq.  (82)  and 
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Flg.    1    Plot  of  log  (K   /T)  vs.    10    /T  for  the  reaction  Br23?:2Br, 
where  K     is  repressed  in  atmos.  ;  data  taken  from  Bureau 
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Fig.   2    Adiabatic  flame  temperature  T^  as  a  function  of  composition 
for  H2-Br2  nnixtures  (p  =  1  atmos.  ,    T     =  323    K) 
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Fig.   3    Equilibrium  mole  fractions  as  a  function  of  gas  composition 

for  H^-Br^  mixtures  (p  =  1  atmos.  ,    T     =  323    K).     The  values 

of  Xtt  c  3-re  so  small  that  they  cannot  be  shown  on  this  plot 
J-l,  I 

(see  Table  III  for  data  on  Xx,  r). 
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48    f— 


Thermal  conductivities  calculated 
without  Eucken  correction  by  Cooley 
and  Anderson  from  the  elementary 
kinetic  theory  of  gases.  (5) 
Thermal  conductivities  calculated 
from  Eq.   (82)  and  the  data  of  Camp- 
bell and  Hirschfelder/^' 
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Fig.   4    The  thermal  conductivity  Ar  as  a  function  of  mixture  connposition 


for  Hy-B^y  fl3-mcSo 


51- 


.5^ 


0 

I— I 

U 


1.0 


1-a 


Fig.    5    Average  specified  heat  c     as  a  function  of  composition  for  the 

ir 


H 


^-Br^  flame  (p  =  1  atmos.  ,    T^  =  323     K).     The  effect  of  die 


sociation  on  c     has  been  neglected. 
P  ^ 
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(5) 
•  Calculated  data  using  thermal  conductivities  of  Cooley  &  Anderson^    ' 

▲  Calculated  data  using  thermal  conductivities  calculated  from  the  data 

of  Campbell  and  Hirschfelder<^) 
^  Range  of  experimental  results  with  burner  size  and  inert  atmosphere 

varied^'*) 
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Fig.   6    Comparison  of  observed  and  calculated  values  for  the  laminar 
burning  velocity  in  H^-Br^  mixtures. 
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